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A comprehensive experimental investigation of exciton and biexciton recombination in 
GaAs/AlGaAs quantum wells is presented. Exciton and biexciton recombination times are found to 
be 16 and 55 ps, respectively. A method of determining the dynamics of the exciton temperature is 
developed. It is shown that both exciton and biexciton recombination processes increase the exciton 
temperature by an amount as high as ~ 10 K. These processes impose a new restriction on the 
possibility of exciton Bose-Einstein condensation and make impossible its achievement in a system 
of direct excitons in GaAs quantum wells even for resonantly excited exciton gas. 
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INTRODUCTION 

Excitons in semiconductors are a simple and very con¬ 
venient analog of an atomic system. Like atoms, exci¬ 
tons can form molecules (biexcitons 0 ) and ions (trions 
0 )- The Rydberg energy and Bohr radius of excitons 
are ~ 10 meV and ~ 10 nm as compared to ~ 10 eV 
and ~ 0.1 nm for atoms. This makes excitons more sus¬ 
ceptible to external fields and allows one to implement 
in the laboratory conditions that are achievable only at 
extreme places of the Universe for atoms. Examples in¬ 
clude the stabilization of the electron-hole liquid by a 
magnetic field Q and spatial separation of electrons and 
holes by an electric field in quantum wells (QWs), which 
makes excitons indirect, suppressing the formation of ex¬ 
citon complexes and enhancing exciton radiative lifetime 


The achievement of Bose-Einstein condensation (BEC) 
in a system of ultracold atoms two decades ago 0, E] to¬ 
gether with rapid progress in nanotechnology inspired the 
work which has led to the attainment of BEC in a system 
of spatially indirect excitons in QWs and exciton 

polaritons in microcavities [12[ [l3| . Atomic BEC was 
achieved at temperatures of ~ 10 - ' K, which comprises 
the main difficulty in the experiment. On the other hand, 
the BEC of excitons can be reached at temperatures of 
~ 1 K due to their much smaller effective mass com¬ 
pared to the atomic mass. Excitons inevitably interact 
with thermal reservoir, the lattice, through the emission 
and absorption of phonons, and exciton system cooling is 
accomplished via the lattice cooling. On the other hand, 
an atomic system can be left on its own in a magneto¬ 
optical trap without interaction with any thermal reser¬ 
voir, and the temperature of such a system is decreased 
by laser and evaporative cooling techniques 14]. Sim¬ 


ilarly to atomic BEC, the achievement of exciton BEC 
suffers from inelastic exciton-exciton collisions leading to 
biexciton formation IT JTb jTol and even to collapse into 
electron-hole liquid [20-221. The exciton system in di¬ 
rect gap semiconductors suffers also from radiative re¬ 
combination. This is why exciton BEC was searched 


for indirect excitons with suppressed recombination and 
Coulomb attraction. 

So far, it was believed that the main effect of exci¬ 
ton radiative recombination and biexciton formation on 
BEC is the decrease of exciton population which can be 
faster than thermalization of the system with the lattice 
0- In this case, it would be possible to attain BEC for 
resonant excitation of cold gas of direct excitons within 
an exciton population lifetime which is ~ 0.5 ns in GaAs 
QWs. This lifetime and even exciton recombination time 
[23I ] are longer than exciton thermalization time 24, 2f| . 
However, so far, only coherence mediated by the exci ting 
laser was observed for exciton gas at low densities 26- 
[32j ], and no signatures of spontaneous coherence were re¬ 
ported for direct excitons, to the best of our knowledge. 
In this paper, we show that exciton and biexciton re¬ 
combination leads to significant heating of the exciton 
system. This imposes a new restriction on the possibil¬ 
ity of exciton BEC. In particular, recombination heating 
makes it impossible to achieve BEC in a system of direct 
excitons in GaAs QWs even for resonantly excited cold 
exciton gas. 

First, we determine the exciton and biexciton recombi¬ 
nation times in QWs. In previous studies, the exciton re¬ 
combination time (radiative decay time of low-momenta 
radiative excitons) was extracted from the photolumi¬ 
nescence (PL) kinetics of resonantly created excitons us¬ 
ing sophisticated models taking into account different 
scattering mechanisms m m. The exciton recombi¬ 
nation time was also determined in four-wave mixing ex¬ 
periments 34J and calculated theoretically 35]. Here, 
by using exciton-resonant excitation with different po¬ 
larizations and powers and also different temperatures, 
we clearly show contributions of different dephasing pro¬ 
cesses to the decay rate of the fast component in the PL 
kinetics and determine the exciton recombination time. 
The biexciton recombination time so far was determined 


only in four-wave mixing experiments [34] and calculated 
theoretically 36, (37]. In this study, we determine biexci¬ 
ton recombination time directly by observing the decay 
of biexciton PL and rise of exciton PL upon the resonant 
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excitation of biexcitons. 

Heating by exciton recombination (evaporative optical 
heating) was considered previously theoretically for free 
carriers [38j and indirect excitons [39] in QWs. How¬ 
ever, to study this effect experimentally, one has to mea¬ 
sure fairly small temperature changes. Here, we develop 
and ground a sensitive method to determine the exciton 
temperature dynamics based on the temperature depen¬ 
dence of the exciton population decay rate. We show 
both theoretically and experimentally that, for the con¬ 
sidered system of direct excitons, the exciton tempera¬ 
ture is higher than the lattice temperature by > 4 K 
due to exciton recombination heating. Furthermore, we 
demonstrate that the formation of biexcitons followed by 
their recombination also contributes significantly to the 
heating of the system. With increasing excitation power, 
the temperature increases by an amount proportional to 
the ratio of the biexciton and exciton concentrations, and 
this amount becomes as high as ~ 10 K. 

EXPERIMENTAL DETAILS 




Energy (eV) t (ps) 

FIG. 1. (a) PL spectra, corresponding to the different times 
after the excitation pulse (symbols) fitted with Lorentzian 
peaks (lines). The spectra are normalized to the maximum 
value and vertically shifted, (b) Dynamics of the exciton (full 
symbols) and biexciton (open symbols) PL intensities. (a),(b) 
Exciton-resonant excitation (A = 0 meV) with P = 0.25 mW, 
Platt = 10 K. 


range. 


The sample under study is a GaAs/Al 0 . 05 Ga 0 . 95 As het¬ 
erostructure with two shallow tunneling-isolated QWs 
of widths d = 3 and 4 nm. The thickness of the 
Alo. 05 Gao. 95 As barrier layer separating the QWs is 
60 nm. Only states of the wider QW were excited and 
studied. The same sample was used in Ref. [4(|, where 
the formation of trions as a result of the preferential cap¬ 
ture of one species of carriers into the QWs under above¬ 
barrier excitation was studied. 

The sample is mounted in a He-vapor optical cryostat 
and excited by the radiation of a mode-locked Ti-sapphire 
laser generating a periodic train of 2.5-ps-long pulses at 
a repetition rate of 76 MHz. The excitation laser beam 
is focused into a 10 - to 20-fj.m spot on the sample sur¬ 
face using a 6 -mm-focus micro-objective located in front 
of the sample surface so that the surface is near its focal 
plane. The PL is collected by the same micro-objective. 
The excitation beam was slightly misaligned with respect 
to the optical axis, which allows to block the direct re¬ 
flection. The micro-objective is mounted on the sample 
holder inside the cryostat, which provides the good sta¬ 
bility of the system against vibrations. The PL coming 
out from the cryostat is focused with a ~ 100 -mm-focus 
lens to form a magnified image of the PL spot on the 
slit of a spectrometer coupled to a Hamamatsu streak 
camera. The slits of the spectrometer and streak camera 
selected the central region of the PL spot with homoge¬ 
neous intensity distribution. In all experiments except 
those described in Section „ (Fig. [2]), the excitation beam 
was linearly polarized and PL was registered in a per¬ 
pendicular linear polarization. The spectral resolution 
is 0.2 — 0.5 meV. The temporal resolution is 5 — 20 ps, 
depending on the spectral resolution and the used time 


RESULTS AND DISCUSSION 


The QW emission spectra for different times t after a 
resonant excitation pulse are presented in Fig. [Ha'). The 
spectra feature two lines separated by about 2 meV. The 
higher-energy line (Ex) is attributed to exciton emission. 
It persists in the spectra even for sufficiently high tem¬ 
peratures and low e-h densities. The lower-energy line 
(BiEx) is significantly reduced with respect to the exci¬ 
ton line as the e-h density decays with time (Fig. [Qa))- 
We attribute this low-energy line to biexciton emission. 
Indeed, generation of biexcitons is the most efficient for 
the resonant excitation [TH], [l(| . The spectrum is fitted 
with two Lorentzian peaks (lines in Fig. [lja)) to deter¬ 
mine the intensities, widths, and spectral positions of the 
exciton and biexciton lines. The intensity of the BiEx 
line decays about two times faster than the intensity of 
the Ex line (Fig. [2(b)) as expected for biexcitons 17. 181] . 
and the biexciton nature of the BiEx line will be further 
confirmed throughout the paper. 


Exciton recombination time 

Here we determine the exciton recombination time (ra¬ 
diative decay time of low-momenta radiative excitons) by 
investigating the exciton dynamics with temporal resolu¬ 
tion of about 5 ps for exciton-resonant linearly polarized 
excitation. Figure [2](a) shows exciton intensity dynamics 
registered in polarization parallel (||, solid curves) and 
perpendicular (_L, dashed curves) to the excitation po¬ 
larization. Two decay components are clearly seen in the 
|| polarization kinetics, while the _L polarization kinetics 
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FIG. 2. Exciton-resonant excitation (A = 0 meV). (a) Exciton intensity dynamics in linear polarizations parallel (solid lines) and 
perpendicular (dashed lines) to the excitation polarization for different excitation powers. TLtt = 10 K. (b) Power dependences 
of the fast component decay time for the linearly (squares) and circularly (circles) polarized excitations. TLtt = 10 K. (c) 
Temperature dependence of the fast component decay time for linearly polarized excitation with P = 0.013 mW (different 
excitation spot position on the sample). 


shows the slow component only. We interpret these ob¬ 
servations as follows. A resonant excitation pulse creates 
a large population of low-fc radiative excitons with wave 
vectors within the light-cone. Accordingly, the fast com¬ 
ponent is induced by the initial low-fc exciton recombi¬ 
nation and scattering towards the nonradiative reservoir, 
while the slow component corresponds to the decay of the 
reservoir exciton population. The contribution of a spin 
relaxation to the decay of the fast component is small 


34]. Otherwise, the PL decay in || polarization would be 


accompanied by the comparable PL rise in the _L polar¬ 
ization, which is not seen in the experiment (Fig. [21(a)). 
The initial scattering of the low-fc excitons to the nonra¬ 
diative states corresponds to the increase of the mean ex¬ 
citon energy (temperature) either due to the absorption 
of phonons (the corresponding scattering rate is l/T p hon) 
or to the formation of biexcitons (the corresponding rate 
is 1/rf). Indeed, each biexciton formation event releases 
energy for the exciton system approximately equal to the 
biexciton binding energy. For the decay rate of the fast 
component neglecting the scattering back from the reser¬ 
voir we can write 1/r = 1/to + l/r p hon + 1/rf. As the 
excitation density P is increased, the rate of biexciton 
formation increases and r should decrease, which we do 
observe in the experiment (black squares in Fig. [21(b)). 
In the limit of P —>■ 0 we have 1/r = l/r 0 + l/rphon- 
Interestingly, for circularly polarized excitation, the fast 
component decay time for the co-polarized intensity is 
almost independent of P. This can be explained by the 
fact that biexciton formation is suppressed for circularly 
polarized exciton population since total biexciton spin 
equals zero fl[ [41 1. 


An increase in the lattice temperature Ti at t should fa¬ 
vor phonon absorption and lead to an increase in l/r p hon, 
causing a decrease in r. However, Fig. [21(c) shows that r 


is almost independent of Ti at t, indicating that r p hon r. 
As we show in the following, interaction with phonons 
is characterized by times of ~ 100 ps. Note also that 
at low P secondary emission can preserve laser-mediated 
coherence for a time comparable to r 0 f2Gi - f32l ] . which, 
of course, does not prevent exciton radiative decay at 
the rate 1/tq. Finally, at small P, r « r 0 and we get 
the exciton recombination time tq ss 16 ± 2 ps. This 
result is consistent with the results of kinetic measure¬ 
ments of Ref. 231 (10 ps for a 4.5-nm QW) and Ref. 


(20 ps for an 8-nm QW), with four-wave mixing stud¬ 
ies of Ref. [34] (13 ps for a 25-nm QW) and theoretical 
calculations Ref. [35( (25 ps for a 10-nrn QW). 


Biexciton recombination time 

Now we determine the biexciton recombination time 
(radiative lifetime) by investigating the exciton-biexciton 
dynamics with a temporal resolution of about 5 ps for 
biexciton-resonant linearly polarized excitation (Fig. [3]). 
Biexciton-resonant excitation implies the absorption of 
two photons leading to the creation of a biexciton hav¬ 
ing an energy of 2E X — Pb- Thus, we used the ex¬ 
citation photon energy Huj = E x — E^/2 detuned by 
A = Piaser — E x = — 1 meV from the exciton resonance. 
Figure [3] shows, that for this excitation energy, the biex¬ 
citon intensity greatly exceeds the exciton intensity at 
short times. Then, the biexciton intensity rapidly decays 
with a characteristic time t\ ss 35 ps. The decay of the 
biexciton intensity is accompanied by an increase in the 
exciton intensity, which starts to decrease at later times. 

The simplest scenario explaining the experimental ob¬ 
servations is the following. An excitation pulse creates 
biexciton population. In a process of biexciton radiative 
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FIG. 3. Dynamics of exciton (squares) and biexciton (cir¬ 
cles) PL intensities for biexciton-resonant excitation (A = — 1 
meV) with P = 0.13 mW. The inset shows power depen¬ 
dences of the exciton (squares) and biexciton (circles) time- 
integrated intensities for biexciton-resonant (full symbols) 
and exciton-resonant (open symbols) excitations. Thtt = 10 
K. 


decay, one of the e-h pairs constituting a biexciton recom¬ 
bines to give a photon with energy E x — E\> and leaves an 
exciton. Thus, the radiative decay of biexciton popula¬ 
tion is accompanied by the creation of the same amount 
of excitons. Subsequently excitons either decay radia- 
tively or form secondary biexcitons giving much weaker 
emission at longer times, so the latter process can be ne¬ 
glected. According to this scenario the time-integrated 
biexciton intensity should be equal to the time-integrated 
exciton intensity. In the experiment, however, the to¬ 
tal exciton intensity exceeds the total biexciton intensity 
by a about factor of r « 3 for all considered excitation 
powers (see the inset in Fig. [3] solid squares and circles, 
respectively). This fact indicates that excitons can be 
formed not only as a result of biexciton radiative decay, 
but also as a result of biexciton dissociation. Further¬ 
more, the excitation pulse can create some initial exci¬ 
ton population which can be seen in Fig. [3] as a nonzero 
value of the exciton intensity corresponding to the maxi¬ 
mal biexciton intensity. The inset in Fig. [3] demonstrates 
that both exciton and biexciton time-integrated intensi¬ 
ties show almost identical power dependences, which at 
low P can be approximated by oc P 119 and oc P 112 , 
respectively, with the exponent about twice as large as 
that for the excitons upon exciton-resonant excitation 
(oc P 0 - 54 , open squares in the inset, the sublinear de¬ 
pendence might be related to the fact that we do not 
register the emission in the direction of the direct reflec¬ 
tion of the resonant beam which shows a large secondary 
emission intensity within a few picoseconds 26J). These 


P-dependences confirm the fact that biexcitons are cre¬ 
ated by two-photon absorption process. Furthermore, 
the possible initial exciton population is also created by 
two-photon rather than by the standard one-photon ab¬ 
sorption process. However, more likely all the excitons 
are created from decay or dissociation of biexcitons only, 
and the initial exciton population reflects the number of 
excitons created from biexcitons within our time resolu¬ 
tion. 

The dynamics of the exciton-biexciton system can be 
described by the following rate equations 17]: 


dN x 

dt 


_N x 

Tx 

dN h 


Nb 

n 


2 ^ ^2 P(iVx), 
Td 


dt =-—- — + F (N*), 

dt Tb T d 


( 1 ) 

( 2 ) 


where N x and N\, are the exciton and biexciton concen¬ 
trations, respectively; r x , Tb, and Td are the times of the 
whole exciton population radiative decay, biexciton re¬ 
combination, and biexciton dissociation, respectively; the 
term F(N X ) describes biexciton formation from two ex¬ 
citons. In the description of the fast stage we disregard 
the formation of secondary biexcitons described by the 
last term in both equations. This process determines the 
biexciton intensity when quasiequilibrium between exci¬ 
ton and biexciton populations is established. This inten¬ 
sity is much smaller than the biexciton intensity at the 
considered fast stage for the same level of the exciton in¬ 
tensity (compare Figs. QJb) and [3J. Equations (HD and 
@ without the last term have the following solution: 


JVx = K ^ ^ [(1 + <*)e-* /T * ~ e-^ 1 ], (3) 

N h = N°e~ tF \ (4) 


where 1/t\ = 1 /t\, +1 /t^ is the biexciton intensity decay 
time, iV 9 is the initial biexciton concentration, parameter 
a > 0 determines the initial exciton concentration: IV 9 = 
aN^(2/ri — 1/Tb)/(1/Ti — 1 /t x ). It is easy to show that 
the ratio of the total time-integrated exciton intensity to 
the total biexciton intensity is 


r = (2— — l)(a:——-b 1), (5) 

n t x - ti 


and, correspondingly, the expression for biexciton recom¬ 
bination time is 


Tb = ^y[l + r(a ^ + 1 ) 4 ]. 

Z T x T\ 


( 6 ) 


Thus, Tb < Ti(l + r)/ 2, and since 1 /ti = 1/Vb + 1/Td, 
Tb > t i. For experimental values t\ = 35 ps and r = 2.9 
we have 35 < Tb < 68 ps. To get the more precise value 
of Tb, we extract a = 0.17 and t x = 69 ps (such a short 
t x is related to the low exciton temperature in the first 
100 ps) from the fit of the exciton intensity kinetics with 
Eq. d3j). Finally we obtain the biexciton recombination 
time Tb = 55 ps and dissociation time Td = 94 ps. 
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Interestingly, the biexciton recombination time Tb = 
55 ps is larger than the exciton recombination time 
t 0 = 16 ps. However, r 0 and Ty have different phys¬ 
ical meanings. The exciton recombination time deter¬ 
mines the rate of radiative decay of only low-fc radiative 
excitons, while the total exciton population, mainly con¬ 
sisting of nonradiative reservoir, decays at longer times 
r x ~ 400 ps, as it follows from the experiment (Fig. |TJb) ) 
at t > 200 ps. On the other hand, biexciton radiative 
decay is allowed for any biexciton wave vector since it 
can be transmitted to the residual exciton. Thus, Tb de¬ 
termines the radiative decay of the whole biexciton pop¬ 
ulation. The value Tb = 55 ps > To is compatible with 
theoretical results of Ref. [36|, where Tb was calculated 
for all biexciton momenta. On the other hand, calcula¬ 
tions performed for biexcitons with almost zero momen¬ 
tum within the giant oscillator strength model 37] as 
well as four-wave mixing experiments, which also probe 
almost zero-momentum biexcitons, [341 ] give somewhat 
smaller biexciton recombination times: 20 and 11 ps, re¬ 
spectively. The obtained value of Tb = 55 ps is much 
smaller than the biexciton recombination time evaluated 
in Ref. [l8| (330 ps) as the time of the biexciton inten¬ 
sity decay under intense nonresonant excitation in the 
assumption that the biexciton is the dominant species. 
As we will show later this assumption is not valid. 


Exciton-biexciton thermodynamics 


Now, let us consider the behavior on the time scale 
t > tq , Tb. One can easily estimate that, contrary to 
the common belief 18], the biexciton population is much 
smaller than the exciton one even for comparable ex¬ 
citon and biexciton intensities I x ~ /b : N d /N x ss 
A)Tb /I x t x ~ Tb/ t x ~ 55/400 ~ 0.1. Furthermore, tak¬ 
ing into account the relatively slow change of the biex¬ 
citon concentration dN^/dt ss 2N\ y /t x at the considered 
time scale (Fig. [l|), it follows from Eqs. (Q]) and (J2J that 
N h » N?/(l + T d /T b ) « 0.41V b eq , where N° q = r d F(N x ) 
is the equilibrium biexciton concentration. Thus, the 
biexciton concentration is far from the equilibrium one 
due to the short radiative decay time of biexcitons com¬ 
pared to the biexciton formation/dissociation time. The 
biexciton population is driven by the excitons just like 
the nonequilibrium polariton population in a microcav¬ 
ity is driven by the exciton reservoir 42]. 

It is instructive to estimate the absolute value of the 
exciton density N x . Poor knowledge of the QW absorp¬ 
tion coefficient makes it difficult to determine N x directly 
from the excitation density. For our estimates, we used 
the exciton line blue-shift which is, according to Ref. [43[ , 
6E ss 3-E x bA| D lV x (the detailed calculations of exciton- 
exciton interaction are given in Refs. [3 SI), where 
E x b is the exciton binding energy and A 2 D is the two- 
dimensional exciton Bohr radius defined as in Ref. [44 1. 


For excitation power P = 0.25 mW (Fig. [T]) at the begin¬ 
ning of the kinetics SE « 0.2 meV; taking E x b = 7 meV, 
A 2 D = 10 nm, we obtain N x « 10 10 cm” 2 . The biexciton 
density can be determined from the exciton density in 
two ways: (i) From the ratio of the biexciton and exciton 
intensities as it was described above: N d ~ N x I d T d /I x T x . 
For N x = 10 10 cm” 2 and Ih/I x = 0.4 at the beginning 
of the kinetics for P = 0.25 mW (Fig. [U we obtain 
N d S3 0.6 x 10 9 cm” 2 , (ii) By calculating the equilibrium 
biexciton density lV b q corresponding to a given exciton 
density from the law of mass action }l8| : 


TV 2 _ g 2 m x T ^ 

? ' req Anh 2 6XP ^h/T), 


K q 


(7) 


where m x S3 0.3 x 10 27 g is the two dimensional exciton 


effective mass [46[, g x = 4 is the exciton spin degeneracy, 
and we put &b = 1. For N x = 10 10 cm” 2 , T = 10 K 
we obtain lV b q w 1.3 x 10 9 cm” 2 . As mentioned above, 
due to the short biexciton lifetime, actual biexciton den¬ 
sity Nb S3 0.4iV b q S3 0.5 x 10 9 cm” 2 . Good agreement 
between biexciton densities calculated in two different 
ways justifies our estimate for the exciton density. 

Although the ratio of the total biexciton and exciton 
concentrations is far from equilibrium, both exciton and 
biexciton momentum distributions should be close to the 
thermal distribution characterized by the internal tem¬ 
perature T in general not equal to the lattice temper¬ 
ature Tj att . This is true since the thermalization time 
determined by interparticle scattering at considered ex¬ 
citation densities is shorter than t 0 and Tb 24, 25j (in¬ 
terparticle scattering time can be estimated from the ho¬ 
mogeneous density-dependent broadening which is up to 
1 meV in our experiments). Short thermalization time 
also excludes the effects related to the laser-mediated co¬ 
herence which may be important only at low densities 
2^-32|. Furthermore, the estimated exciton density in¬ 
dicates that the exciton ensemble is nondegenerate at 
the considered temperatures and can be described by the 
classical Boltzmann distribution. 


Internal temperature dynamics 

Now, we will describe how the internal temperature T 
of the exciton-biexciton system can be determined from 
the decay time of the exciton population t x , and discuss 
the dynamics of T. 

The radiative decay rate of the whole exciton pop¬ 
ulation is given by the following relation: N x /t x = 
Aj-ad /to, where A rad is the concentration of excitons in 
the radiative zone with wave vectors within the light 
cone k < ojn/c and spins ±1. This concentration can 
be evaluated as A rad ss SE r (l/2)f(0)D = SE r N x /2T , 
where we put the energy of the bottom of the exci¬ 
ton dispersion curve to zero, SE r = (frion / c) 2 / 2m x ~ 
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FIG. 4. (a) Dynamics of the exciton PL intensities (symbols) 
at different lattice temperatures for exciton-resonant excita¬ 
tion with P = 0.13 mW. The lines show the exponential fits, 
(b) Lattice temperature dependence of the exciton PL decay 
time (symbols). The line shows the linear fit. 


80 /xeV is the maximal energy of the radiative excitons, 
f(E) = (Trh 2 N x /2m x T)exp(—E/T) is the occupancy 
of exciton states given by the Boltzmann distribution, 
D = 2m x /Trh 2 is the exciton density of states per unit 
area, n ~ 3.3 is the refractive index of the GaAs, and we 
used T » SE r . Finally, for the radiative decay time of 
the whole exciton population we have 

r x = 2t 0 T/5E, = j, (8) 

and we introduced phenomenological constant A = 
6E t /2tq which will be determined in the experiment and 
related to the measured To. 

Figure 0 a) shows the exciton PL kinetics at different 
lattice temperatures. At long times, the PL decay can be 
described by a monoexponential function and the decay 
time is equal r x if we neglect biexciton formation. The 
dependence of r x on Ti at t is shown in Fig. 0b). As ex¬ 
pected, this dependence is linear r x = (Ti at t + T lec )/A , 
however, there is nonzero offset suggesting that the in¬ 
ternal exciton temperature is higher than lattice temper¬ 


ature by a constant value T rec ss 4 K. The origin of T rec 
will be discussed in the following. The experimental value 
A = 37 K/ns corresponds to the exciton recombination 
time To « 13 ps in good agreement with measured value 
of 16 ps. The linear dependence of the exciton population 
decay time on temperature was also observed in previous 
experiments [Ref. [47[ (A ss 20 K/ns for a 4-nm QW) and 
Ref. 0 {A ss 45 K/ns for a 4-nm QW)] and calculated 
theoretically [Ref. [35] (A « 29 K/ns for a 10-ran QW)]. 

We have shown that the exciton internal temperature 
is proportional to the decay time of the exciton intensity 
(Eq. ®). However, in the experiment the decay time can 
be defined unambiguously only at large times, when the 
intensity decay is exponential and the internal tempera¬ 
ture almost reached its steady-state value. Now, we ex¬ 
tract the internal temperature at any time T(t ) from the 
dynamics of the exciton and biexciton intensities. The 
exciton intensity (defined as the number of particles emit¬ 
ted in the unit time) I x (t) = N x (t)/r x (t) = AN x {t)/T{t). 
On the other hand, N x (t) = / ( °° \l x (t') + Ib{t'))dt' , where 
we use N\, <C N x . Thus, 


T(t) = A 


Ut) 


(9) 


This expression offers a method to determine the exciton 
temperature from the PL dynamics. In the case of mo¬ 
noexponential decay and when Ib(t) <C / x (t), Eq. (0 is 
reduced to Eq. ©• 

To test the validity of this method, we perform ex¬ 
periments using excitations with different excess ener¬ 
gies A = £\ aser — E x with respect to the exciton res¬ 
onance, which defines the initial exciton temperature. 
The excitation density is relatively low, so the biexci¬ 
ton recombination heating (see following) is small (for 
A < 0 it is small at t > Tb). The exciton PL dynamics 
under excitation with different excess energies is shown 
in Fig. 0a). Figure 0b) shows the exciton tempera¬ 
ture dynamics determined according to Eq. 0 from the 
measured exciton and biexciton intensity dynamics. As 
expected, the temperature relaxes from the value deter¬ 
mined by the excess energy and reaches the steady-state 
value T st = TJ att + T rec at long times. The following qual¬ 
itative criterion applies: the faster the intensity decrease, 
the lower the temperature and vice versa [compare curves 
in Figs. 0a) and0b)]. Our method can be applied even 
for low exciton temperatures and concentrations, where 
the temperature can not be determined from the Boltz¬ 
mann tail in the PL spectra, which originates either from 
the e-li plasma recombination 0] or from the longitudi¬ 
nal optical (LO) phonon replica of exciton recombination 
0 . 


Now, we show that increase of the internal exciton tem¬ 
perature T by T rec [offset of the de pen dence in Fig. 0b)] 
is related to recombination heating [38( (evaporative opti¬ 
cal heating Q), a process similar to evaporative cooling 
in atomic systems [8|, [14]. The average exciton energy is 
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FIG. 5. Dynamics of the exciton intensity (a) and internal 
temperature (b) for excitations with different excess energies 
A with respect to the exciton resonance and P = 0.13 mW. 
Lines in Fig. [5jb) show exponential fits. Ti at t = 10 K. 


T, but only excitons with low energies < SE r -C T recom¬ 
bine; as a result, the average exciton energy is increased. 
Here, we disregard biexciton formation, which is insignif¬ 
icant at long times. To describe exciton recombination 
heating quantitatively, let us consider the evolution of 
the exciton system total energy per unit area N X T: 


from rising infinitely. By fitting the temperature dynam¬ 
ics in Fig.[5jb) with Eq. (fill) . we determine k « 5.4 ns 1 
(exciton-phonon scattering time 1/k « 180 ps), and 
we can calculate the recombination heating temperature 
T rec = A/k w 37/5.4 ss 6.9 K which is in a reasonable 
agreement with the value T rec = 4 K determined from 
the data in Fig. Q|b). 

Another process that, in addition to exciton recombi¬ 
nation heating, contributes to the increase in the exciton 
temperature is the formation of biexcitons with their sub¬ 
sequent recombination. Indeed, each biexciton formation 
event with subsequent recombination releases energy E b 
for the exciton system, while the emitted photon has en¬ 
ergy — E b . Thus, taking into account that N h <C AT X , 
while intensities 4> = N b/ib and 4 = N x /t x might be 
comparable, one has to add a term EbNb/rb to the right- 
hand side of Eq. to- This term describes the energy 
increase due to biexciton recombination. As a result, the 
temperature dynamics is described by the equation 


dT 

dt 


— A — k(T — Tjatt) + 


Tb Nb 
Tb Ax ’ 


(14) 


leading to 


±(N x T) = -^^-k(T-T um )N x . ( 10 ) 

dt r x 2 

The total energy changes due to the radiative decay (first 
term, where E T /2 is the average energy of recombining 
excitons) and due to the emission/absorption of phonons 
(second term). The latter process tends to equalize the 
temperatures of the exciton system and the lattice, so in 
the first-order approximation, its rate is proportional to 
T — Tiatt- The proportionality coefficient n characterizes 
the strength of the exciton-phonon interaction. Taking 
into account that dN x /dt = —N x /t x = —AN X /T and 
SE r -C T, we obtain the solution of Eq. TO: 

T = T 0 e~ Kt + (Tiatt + -)(l-e- Kt ), (11) 

K 

where To is the initial temperature. Thus, the steady- 
state value of the exciton temperature achieved at t ^> 
1/k 


T s t — T^tt H- 

K 


( 12 ) 


is higher than Ti a tt by a constant T rec = A/k, in agree¬ 
ment with the results in Figs.[4jb) and[5](b). The recom¬ 
bination heating temperature 


Tree - - 


SE r 

2tqk 


(13) 


is determined by the exciton recombination rate and 
exciton-phonon interaction constant. Interestingly, only 
interaction with the lattice prevents the temperature 


T = T 0 e~ Kt + (T^tt + -)(1 - e~ Kt ) 


+ 


K 

Eb 

Tb 


l (15) 


The last term determines the temperature increase due 
to biexciton recombination. To calculate this term, one 
should obtain the time dependence of N x and Nb from 
Eqs. ([J) and ([2]). However, the constants r x , Td and for¬ 
mation rate F(N X ) in these equations in turn depend on 
T. To get an approximate expression for the last term in 
Eq. (fT5]> - we P ut A/x oc exp(— t/r x ) and Nb oc exp(— 2 t/r x ) 
and disregard the dependence of r x on T. Finally, 


T = T 0 e~ Kt + (Tiatt + -)(1 - e ~ Kt ) 

K 


+ - £t , (“> 

TbK - T b /r x / x (0) 


Thus, biexciton recombination can increase the internal 
temperature by as much as ~ Eb- 

Figure [6jb) shows the exciton temperature dynam¬ 
ics for different powers of exciton-resonant excitation. 
The corresponding exciton and biexciton intensities are 
shown in Fig. Ha). For the lowest excitation power, the 
temperature monotonously increases to the steady-state 
value T st . On the other hand, for higher powers, when 
the biexciton intensity becomes comparable to the exci¬ 
ton intensity (Fig. H a ))> the temperature first increases 
above T st and then decreases to T st . This initial temper¬ 
ature rise results from biexciton recombination heating 
and is proportional to 4/4- As the biexciton intensity 








FIG. 6. Dynamics of the exciton (solid symbols) and biexci- 
ton (open symbols) intensities (a) and the internal tempera¬ 
ture (b) for exciton-resonant excitation (A = 0) with different 
powers. Lines in Fig. [6jb) show biexponential fits. Tiatt = 10 

K. 


becomes much lower than the exciton intensity, the tem¬ 
perature relaxes to T s t. The dependences in Fig. [61(b) are 
fitted with double-exponential functions in accordance 
with Eq. (usd with the temperature decay time set to 
r x = 0.38 ns. 

The described heating mechanisms cancel out the pos¬ 
sibility of exciton BEC in the considered system even for 
resonant excitation. Indeed, the system temperature can 
not be lower than ~ 4 K for a time longer than few tens 
of picoseconds even for Ti att = 0. This already gives a 
BEC threshold density of at least _/V t hr ~ 2m x T Tec /n ~ 
1.6 x 10 12 cm~ 2 which is of the same level as the exciton 
saturation density N sat ~ 1 /ag ~ 10 12 cm -2 . Further¬ 
more, an increase in the exciton density leads to a pro¬ 
portional increase in the temperature due to biexciton 
recombination heating Tf, oc N\,/N x oc N x . Experiment 
shows that even a moderate increase in N x , to the level 
when 7b becomes comparable to 7 X , leads to an increase 
in T up to ~ 10 K. This effect forbids exciton BEC even 
if one forgets about exciton saturation. 


the thermodynamics of the whole system is that they 
provide an additional radiative channel for excitons. 

A method of determining the temperature dynamics of 
the exciton system has been developed. It is based on the 
proportionality of the radiative decay time of the whole 
exciton population and temperature of the exciton sys¬ 
tem and also takes into account biexciton recombination 
channel. The method has been tested in experiments 
with different excess energies of the excitation photons 
which define the initial temperatures of the system. 

The heating of the exciton system by exciton recom¬ 
bination has been analyzed theoretically and revealed in 
the experiment. This effect is analogous to evaporative 
cooling in atomic systems. An increase in the tempera¬ 
ture due to recombination heating (~ 4 K in the QWs 
under study) is proportional to the exciton recombination 
rate and inversely proportional to the rate of cooling by 
acoustic phonons. 

If has been shown that biexciton formation and sub¬ 
sequent recombination also lead to an increase in the 
exciton temperature by an amount proportional to the 
ratio of the biexciton and exciton emission intensities. 
For comparable biexciton and exciton intensities, this in¬ 
crease can be as large as tens of Kelvins. 

The revealed recombination heating effects imply new 
restrictions on the Bose-Einstein condensation of exci¬ 
tons. These restrictions cancel out BEC in many ex¬ 
citon systems even where the thermalization rate ex¬ 
ceeds the radiative decay rate (like direct excitons in 
GaAs/AlGaAs QWs) and even for resonantly created 
cold exciton gas. 

The main quantitative conclusions of the paper are 
summarized in Table [U 

We are grateful to M. L. Skorikov for careful reading 
of the manuscript, valuable remarks and advices, and to 
N. N. Sibeldin and V. A. Tsvetkov for useful discussions. 
The work is supported by the Russian Science Founda¬ 
tion (Project No. 14-12-01425). 


CONCLUSION 

A consistent picture of exciton and biexciton recom¬ 
bination in GaAs/AlGaAs quantum wells has been de¬ 
veloped. The exciton recombination time has been de¬ 
termined to be To = 16 ps, while the radiative decay 
time of the whole exciton population including large- 
momentum nonradiative part is r x = T/(37 K/ns) tq- 
The biexciton recombination time, which also character¬ 
izes the radiative decay of the whole biexciton popula¬ 
tion, has been determined to be Tb = 55 ps. As a result 
of Tb <C t x , the biexciton concentration is significantly 
lower than one would expect in thermal equilibrium. It 
is also much lower than the exciton concentration, while 
the exciton and biexciton emission intensities can be com¬ 
parable. Thus, the main role of the biexciton states in 
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Exciton recombination time 

to = 

16 ps 


Biexciton recombination time 

Tb = 

55 ps 


Radiative decay time of the whole exciton population 

Tx = 

T/A, A 

= 37 K/ns 

Exciton temperature dynamics 

T(t) 

=A[/r 

(W) + I h {t'))dt']/Ut) 

Exciton recombination heating temperature 

Tree 

= A/k r 

i 4 K 

Biexciton recombination heating temperature 

Tb = 

: (Eh/rb) Jn[Nb{t')/N x (t')] exp[—K,(t - t')]dt' 


TABLE I. Summary of the main quantitative results. 
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